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ABSTRACT 

For many years NASA has been interested in the storage and transfer of cryogenic fuels in space. 
Lunar, L2 and other chemical propulsive space vehicle missions now have staged refueling needs 
that a fuel depot would satisfy. The depot considered is located in lower earth orbit. Many 
considerations must go into designing and building such a station. Multi-layer insulation 
systems, thermal shielding and low conductive structural supports are the principal means of 
protecting the system from excessive heat loss due to boiloff. 

This study focuses on the thermal losses associated with storing LH2 in a passively cooled fuel 
depot in a lower earth equatorial orbit. The corresponding examination looks at several 
configurations of the fuel depot. An analytical model has been developed to determine the 
thermal advantages and disadvantages of three different fuel depot configurations. Each of the 
systems consists of three Boeing rocket bodies arranged in various configurations. The first two 
configurations are gravity gradient stabilized while the third one is a spin-stabilized concept. 

Each concept was chosen for self-righting capabilities as well as the fuel settling capabilities, 
however the purpose of this paper is to prove which of the three concepts is the most efficient 
passively cooled system. 

The specific areas to be discussed are the heating time from the fusion temperature to the 
vaporization temperature and the amount of boiloff for a specific number of orbits. Each of the 
previous points is compared using various sun exposed surface areas of the tanks. 


INTRODUCTION 

For many years NASA has been interested in the storage and transfer of cryogenic fuels in space. 
Lunar, L2 and other chemical propulsive space vehicle missions now have staged refueling needs 
that a fuel depot would satisfy. The depot considered is located in lower earth orbit. Many 
considerations must go into designing and building such a station. Multi-layer insulation 



systems, thermal shielding and low conductive structural supports are the pnncipal means of 
protecting the system from excessive heat loss due to boiloff 

This study will focus on the thermal losses associated with storing liquid hydrogen (LH 2 ) in a 
passively cooled fuel depot in a lower earth equatorial orbit. Three candidate designs will be 
analyzed and discussed, each resulting from a different depot configuration. Each design is 
based on a Boeing Delta IV Heavy Rocket. The specific analyses used to evaluate the designs 
are the heating time from the fusion temperature to the vaporization temperature and the boiloff 
losses. Multi-layer insulation systems, thermal shielding and low conductive structural supports 
are the principal passive means of reducing boiloff in systems storing cryogenic fuels. 

An extensive analytical model has been developed and tested to determine the thermal losses for 
the three different fuel depot configurations. Each system is shown in figure 1 
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Figure 1 : Abacus Tail Feather Concept, b) Four Wing Concept and c) Spinner Concept 


ANALYTICAL MODEL 

The vaporization temperature for LH 2 is 20. 4K; the corresponding fusion temperature is 14K. 
Storage of LH 2 within this range will prevent slushing and boiloff. The greater concern is 
boiloff. Boiloff is mass loss due to phase changes from a solid to a gaseous state. Boiloff causes 
two severe problems within the system: (1) fuel loss and (2) over pressurization of the tank. Fuel 
loss due to vaporization is very expensive considering the mass transfer rate to LEO is 
$ 10,000/lb. Over pressurization of the tank from fuel vaporization is due density change in the 
fuel. This increase in volume could rupture the tank. 


Analysis Assumptions 

Design of the model required the development of several assumptions to govern the process: 
They are as follows: 

1. All potential and kinetic energy changes are considered negligible 

2. Neglect the heating of the boiloff vapor, it is small enough too become negligible. 

3. Storage depot should maintain an Lower earth equatorial orbit, 

4. Assume steady state surface conditions of the following 

(a) 1/10 of tank is considered against dark space 
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(b) 2/5 of tank is covered by the Earth 

(c) The Sun and the thermal shield are combined to cover !4 the area of the tank 
5. Isothermal conditions for each area section of the tank 

CONTROL VOLUME 

Analysis of the storage tank shown in Eigure 2 begins by the application of the first law of 
thermodynamics to the control volume shown. The control volume in this case is the cryogenic 
tank itself as shown below. 
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Figure 2: Control volume of the system 


This is a passive system, meaning that no additional heat comes in or is taken away by thermal 
control system. The corresponding energy balance equation is E in dot + E g dot + E out dot = E st dot 
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Figure 3: Theoretical control volume of the system 


where E m is the energy into the system from outside sources, E out is the energy loss from the 
system. Each of these energy components is associated only with the surface of the control 
volume, meaning the radiation added and lost to the system. E g is the thermal energy generation 

and E st is the rate of energy stored within the system. The energy stored is in the form of 
increasing LH 2 temperature. 

The control volume and energy balance theory stated above is attributed to the specific problem 
of the cryogenic fuel tank. A orbiting space vehicle experiences three sources of heat 
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addition, E in , radiation from the sun, radiation from the earth and the radiation from the vehicle 
itself The vehicle radiation in this case is from the thermal shield and or support structure. 

The core problem of cooling cryogenic fluids is to combat the radiant heating by direct and 
indirect means. The effective radiation heat transfer constant used for the sun is 1143 W/m 2 
Radiation emitted from the earth at LEO of 400Km is given by the radiation constant of q ear th = 60 
W/m 2 for the lighted portion and q ea rth = 20 W/m 2 for the shadow. Radiation absorption occurs 
from the thermal shield also. The thermal shield blocks some % of the sunlight from direct heat 
transfer to the tanks but in turn the shield is heated and thus emits radiation. The radiation 
constant is assumed to be 10 W/m 2 . 

E m results from the sum of the sources: 


I' in QradSun Q radEarth Qi 


'radShield 


( 1 ) 


Heat is also lost from the system from the tank surface. The effective radiation from the tank 
takes the form as radiation: 


QmdTank ~ ^ ' O’ • A • 7i, 


( 2 ) 


Where e is the emmissivity and a is the Stefan-Boltzman constant. Other heat transfer 
mechanisms and heat lost due to the heat sink in boiloff are neglected: 

R out Qradfank (^) 

Figure 4 depicts the control volume of the tank with the heat additions and losses due to radiation 
heat transfer: 

The total heat balance for the surface of the control volume is: 
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radTot QradSun Q radEarth 
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( 4 ) 


SURFACE TEMPERATURE CALCULATION 

To calculate the surface temperature of the tanks requires the determination of the roots of a forth 
degree polynomial. Equating the conduction and radiation heat transfer rates and the polynomial 
is shown in Eqn 9-12. 
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Figure 4: Control volume for radiation absorbed by the tank from the sun, earth and the shield and 
radiation emitted by the tank 


The properties of the tank skin and substructure affect the thermal heating and cooling. The 
assumed tank absorbivity and emmissivity are assumed a = 0.4 and e - 0.018. Another strong 
factor contributing to the heating and cooling of the cryogenic tanks is the time spent in the sun. 
The assumptions are made that the time spent in earth’s shadow is equivalent to 40% of the time 
in orbit, as shown below. 


TANK CONDUCTION HEAT TRANSFER 


Now that the surface heat exchange has been established the thermal resistance throught the tank 
wall must be established. Thermal resistance is the ratio of the driving potential to the 
corresponding transfer rate. The driving potential for heat transfer through the tank is the outside 
surface temperature inside fuel temperature, (Tj -T s ). Thermal resistance for conduction through 
multi-layered insulation (MLI) is modeled as follows: 


^Fuel — 


Rj - ■ La 


,R„=- Lp 


huefA’“ A k A -A ’ F k P - A 


(9-11) 


Where h is the convection heat transfer coefficient, k is the thermal conductivity and L is the wall 
thickness. 


Equivalent thermal resistance through tank wall: 


R-EQ — R ! ‘'ud + R/i + R/‘ 


( 12 ) 


The resulting conductance equation is as follows: 
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Cond 


Req 


(13) 


A visual model of the tank wall and the corresponding resistance is seen below. 
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Figure 5: Equivalent thermal circuit for a series resistance wall of the Delta IV Tank 


BOILOFF 


The description for the rate at which energy is stored in the system is in terms of the temperature 
increase in the fuel and tank wall. This phenomenon occurs succeeding the radiation transfer due 
to small transport delays. Once the inner surface of the tank, T s , changes the results are directly 
related to the changing temperature in the fuel, through natural convection heat transfer. The 
nucleate pool boiling equation to evaluate the heat transfer to the fuel from the tank wall is as 
follows: 


9, 
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(14) 


Where A is area, C s fis the coefficient for surface tension combination, Cpis the specific heat at 
constant pressure, g gravity at lower earth orbit, hf g latent heat of vaporization, P r ; s the Prandtl 
number, p is theViscosity, p is density of liquid fuel and a is the Stefan-Boltzman constant. 

Once the boiling heat transfer rate has been determined, assuming temperature and surface area 
does not change, all heat addition to the tanks will result in raising the temperature and or boiloff 
Determination of the boiloff rate is accomplished using Eq. 15. 


m 


b,dot 



q s is the boiling heat transfer rate and latent heat of vaporization. 


(15) 


MODELS 

As mentioned this paper addresses three different types of models, abacus tail feather concept, 
four wing concept and spinner concept. Each of the three concepts possesses certain similar 
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thermal characteristics. All three of the concepts are based on construction of Delta IV Boeing 
rocket tanks. Thus, the thermal properties and geometry of the tanks are the same. In addition to 
the thermal and geometric similarities of the concepts, a thermal shield is used with each design. 

The shield and the sun combined are assumed to cover approximately 50% of the total area of the 
fuel tanks, while the remaining 50% of the surface area is exposed to either the earth or space. 
The ratio of sun-exposed area to the shaded area depends solely on the area of the tank in direct 
sunlight. For example, if the sun hits 10% of the total area, the shield will shade the remaining 
40%. This is true for the abacus and four wing concepts. An exception, however is the spinning 
concept, for which the amount of sunlight is dependent on rotation period of the tanks. 

The abacus concept has three tanks situated at the bottom of the thermal shield. The shield will 
cover each tank for a portion of the time in orbit. Forty-percent of the orbit time is spent in the 
shadow of the earth and the remaining sixty-percent of the time is a constant exposure to 
sunlight. The percentage of exposure can vary from 2% to 50% of the total area of the tanks. 
Naturally, the less time in the sun the lower the heat gains. 

The spinner concept also consists of three tanks. The tanks are attached to a central hub located 
beneath the thermal shield. The tanks are spun for cooling and fluid collation purposes. The 
spin rate of the tanks to this point is not valued; however, an assumed value is used here for 
calculation purposes. This value is calculated as follows 

Acceleration due to spinning is: 

a = co 2 ■ r = coxcox(-r) (16) 

Solving for get an effective artificial gravity of % g gives: 

co = b rad / sec (1 7) 

Thus the total rotation speed for the system is: 

T ro,atio n al=Y^ Hz ( 18 ) 

CO 

The tail feather / gravity gradient concept uses a series of arrays mounted on the system. Even if 
the solar cells double as a thermal shield, it will not likely cover a large enough area to 
effectively block the sun’s radiation heat. 


RESULTS 

The gravity gradient is surprisingly the design with the most massive thermal penalties. The 
spinner concept proved to be the best design for reducing mass boiloff in a passively cooled 
system. The abacus and spinner concepts were compared directly using the same time in 
equatorial orbit, and altitude. In addition, the concepts varied exposure to the sun from 2% to 
50% of the total surface area. However, the spinner concept had an opportunity to cool the tanks 
during the penod of rotation in the shade while the abacus ant gravity gradient did not. 
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From figure 5 the heating time from the fusion temperature to vaporization temperature is shown 
for the three concepts for different sun exposed surface areas. The spinner concept is the best for 
the majority of the time, however the abacus method the better for the smaller exposed surface 
area. And the tail feather concept is better for the greatest exposed surface area. 


Depot Heating Comparison 

Abacus » Spinner Four Wing 



Figure 6: Fuel depot heating comparison chart. At any given sun exposed surface percentage the 
heating time from fusion temperature to vaporization temperature is known. 

The analysis of boiloff rates for the three designs showed some staggering results. The spinner 
concept is by far the best design thermally of the three proposed, followed by the abacus concept 
and, last, the gravity gradient concept. The abacus concept had a boiloff amount of 1 19 times 
greater than the spinner concept for 450 orbits. From this analysis it is apparent that a system of 
these types cannot exist with passive cooling in orbit, active cooling must be employed. 


CONCLUSION 

Consideration must be taken for the energy balance using the first law of thermodynamics, hi the 
design of cryogenic storage the minimum boiloff rates of systems in LEO should be optimized. 
Several recommendations for making this project work effectively are: 

1 . Increase the insulation thickness on the Delta IV tanks. Do not shave any off for 
weight considerations; allow the insulation system to be the payload itself, 

2. Decrease the time the fuel depot spends in the sun while subsequently increasing 
the thermally shielded surface 

3. Attach an active cooling system to reduce the effective boiloff 
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4. Use the boiloff gasses to cool the thermal shield, cooling the tank surface and 
possibly, use to cool the MLI 

5. Use a spherical storage; surface tension dictates that a body of fluid neutrally 
floating in a gas will assume the smallest shape possible. Thus a spherical tank 
has the smallest effective surface area. 

6. Design the system such that it has low conductive structural supports. 
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NOMENCLATURE, ACRONYMS, ABBREVIATIONS (STYLE=PAPERHEADINGS) 

A area, m 2 

a absorbivity of the fuel tanks 

C s f coefficient for surface tension combination 

Cp specific heat at constant pressure, J/kg*K 

e emissivity of the fuel tanks W/m 2 *K 

g gravity at lower earth orbit, m/s 2 

hfuei convection heat transfer coefficient 

hf g latent heat of vaporization, J/kg 

kA thermal conductivity, W/m*K 

L f thickness ofinsulation, m 

mdot mass boiloff rate, kg/sec 

n exponent for surface tension combination 

Orbit number of orbits in analysis 

P r Prandtl number 

Q overall heat delivery to the system, W 

q Radiation constant, W/m 2 

R Resistance 

T Temperature, K 

p Viscosity, kg/s*m 

pi density of liquid fuel, kg/m'' 

p g density of gaseous fuel, kg/m' 

a Stefan-Boltzman constant,W/m 2 *K 4 
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